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Abstract 

The interaction between magnetic impurities and the gapless surface state is of 
critical importance for realizing novel quantum phenomena and new functionalities 
in topological insulators. By combining angle-resolved photoemission spectroscopic 
experiments with density functional theory calculations, we show that surface deposition 
of Cr atoms on Bi 2 Se 3 does not lead to gap opening of the surface state at the Dirac 
point, indicating the absence of long-range out-of-plane ferromagnetism down to our 
measurement temperature of 15 K. This is in sharp contrast to bulk Cr doping, and 
the origin is attributed to different Cr occupation sites. These results highlight the 
importance of nanoscale configuration of doped magnetic impurities in determining the 
electronic and magnetic properties of topological insulators. 

Three dimensional (3D) topological insulator (TI), a novel class of materials,® 2 is a 
topologically nontrivial bulk insulator with conducting surface state (SS). Protected by the 
time-reversal symmetry (TRS), the SS is spin-polarized, gapless and robust against nonmag¬ 
netic impurities.^'" 6 Meanwhile, magnetic impurities may induce out-of-plane ferromagnetism 
and break the TRS in TI, resulting in a band gap at the Dirac point of the SS. 7 Such 
TRS breaking is essential for the realization of various novel quantum phenomena which 
have great application potential in future nano-electronic and spintronic devices, such as 
quantum anomalous Hall effect (QAHE),®^ topological magnetoelectric effect (TME) 12 and 
birefringent spin lens. — Thus, it is critical to understand the interaction between 3D TIs and 
magnetic impurities both for fundamental physics and nano-device applications. 

The interaction between 3D TI and doped magnetic impurities can modify the electronic 
and magnetic properties of 3D TI, and the results may strongly depend on the types of 
magnetic atoms, occupation sites of the magnetic impurities 14 and experimental conditions 
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under which they are introduced. For example, although doping Fe and Cr into the bulk 
Bi 2 Se 3 crystal or thin film during the growth process leads to a gap opening at the Dirac point 
suggesting TRS breaking,^ EHI3 the robustness of the gapless SS against Fe deposition on the 
surface has been debated® 19 and surface deposition of Cr still remains to be investigated. 
Here we present a combined experimental and theoretical study of the electronic structures 
of Bi 2 Se 3 upon surface deposition of magnetic Cr atoms on freshly cleaved B^Ses single 
crystals. Angle-resolved photoemission spectroscopy (ARPES) data show that the SS is robust 
upon surface deposition of Cr, and the absence of gap opening at the Dirac point suggests 
that there is no long-range out-of-plane ferromagnetism down to the lowest measurement 
temperature of 15 K. In addition, Cr atoms dope electrons to the bulk states of BDSes 
forming a two-dimensional electron gas with large Rashba-splitting on the surface. These 
observations are in striking contrast to the bulk doping of Cr during thin film growth, 15 
suggesting that Cr atoms deposited on the surface interact very differently with the TI 
sample. Density functional theory (DFT) calculations show that the difference in electronic 
and magnetic properties originates from the different Cr configurations: hollow site for surface 
deposition versus substitutional site for bulk doping. Due to the coupling between Cr 3d 
orbitals around the Fermi level, the magnetic easy axis of deposited Cr atoms is likely along 
the in-plane direction, which is in agreement with the robustness of the SS. These results 
highlight the importance of nanoscale configuration of magnetic impurities in determining 
the electronic and magnetic properties of topological insulators, and such understanding is 
critical for manipulating the electronic properties of TI and future electronic applications. 

Figure [l] shows a comparison of ARPES data measured at 15 K with 35 eV photon energy 
on freshly cleaved BDSesfTll) surface (a) and after deposition of Cr atoms on the surface for 
46 minutes, with coverage ~ 9% monoatomic layer (ML) Cr (d). ARPES data on pristine 
Bi 2 Se 3 [Fig. m shows conical dispersion from the SS with Dirac point energy at -0.3 eV, a 
small parabolic bulk conduction band near E^. After deposition of Cr for 46 minutes, the SS 
persists and is electron doped, with the Dirac point energy shifted down to -0.51 eV. Energy 
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Figure 1: (Color online) (a) Dispersion image along M-F-M direction. Dotted lines show 
dispersions extracted from MDCs in (c). (b) EDCs at momentum positions between ki and 
k 2 i marked by ticks in (a). The EDC at the F point is highlighted by red color, (c) MDCs 
at energy positions marked by ticks between Ei and E 21 in (a). Ticks are guides for peak 
positions of the SS. (d) Dispersion image along M-T-M direction after surface deposition of 
Cr for 46 minutes, ~ 9% ML Cr on Bi 2 Se 3 . Dotted lines show dispersions extracted from 
MDCs in (f). (e) EDCs at momentum positions between ki and k 29 marked by ticks in (d). 
The EDC at T point is highlighted by red color, (f) MDCs at energy positions marked by 
ticks between Ei and E 24 in (d). Ticks are guides for peak positions of the SS and Rashba 
splitting bands. 
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distribution curves (EDCs) and momentum distribution curves (MDCs) analysis shows that 
there is no gap opening at the Dirac point of the SS down to 15 K [Fig. |TJ[e, f)]. This is in 
striking difference from bulk doping of Cr into B^Ses Elms during MBE growth, 15 20 where 
a significant gap at the Dirac point was reported at a much higher temperature. Despite the 
magnetic moment of Cr atoms®®] ^ our results suggest the surface state is very robust upon 
Cr surface deposition. 




Deposition time (mins) 




Deposition time (mins) 


Figure 2: (Color online) (a-j) Dispersions along M-T-M direction on pristine Bi 2 Se 3 (a) and 
with different Cr deposition time 5, 10, 15, 20, 25, 32, 39, 46 and 57 minutes (b-j). Red 
dotted lines indicate the energy positions of the Dirac point, (k) Evolution of the Dirac point 
energy as a function of Cr deposition time. The inset shows a schematic illustration of the 
band structure evolution. (1) Rashba splitting A kp (black symbols, right axis) and calculated 
Rashba coupling parameter a (blue symbols, left axis) as a function of Cr deposition time. 


While doping electrons to the SS with preserved TRS, Cr deposition also induces additional 
bands (pointed to by arrows in Fig. 1(d)), which are attributed to two dimensional electron 
gas (2DEG) formed on the sample surface. Figures |2]( a-j) show the evolution of the band 
dispersion with increasing Cr deposition time. Initial deposition induces electron doping 
to the sample and shifts the SS and bulk bands to higher binding energy. A parabolic 
band of 2DEG forms at the bottom of the bulk conduction band after depositing Cr for 10 
minutes [Fig. |2j[c)]. With increasing Cr deposition, the 2DEG band becomes more clear and 
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starts to show momentum splitting [Fig. |2][e-j)] due to the formation of effective potential 
gradient across the sample surface. Such momentum-displaced Rashba effect derives from the 
spin-orbit coupling due to the effective potential gradient on the surface, (k x z) ■ cr, where & 

is the Pauli matrices. The dispersion of the Rashba splitting bands can be approximated 

n 2 k 2 

as E ± (ku) = Eq + 2^1 ± a\k\\ |, where m* is the effective mass and a is the Rashba coupling 
parameter. Similar effect has been observed for various gas absorption or metal deposition on 
Bi 2 Se 3 surface, e.g. CO 24 and water, 25 nonmagnetic alkaline metals. 26 From the evolution 
of electronic structure with Cr deposition, we extract the energy position of Dirac point 
E_d [Fig. [2][k)] and the momentum splitting value A kp [Fig. [2](1)] after different deposition 
time. Comparing the evolution of Ejj and A kp, we can see the Dirac point shifts downward 
quickly before the 2DEG band starts splitting after around 20 minutes of Cr deposition. This 
implies the deposited Cr atoms first donate electrons to B^Sea sample and then form 2DEG 
with increasing effective electrical potential near the sample surface. The amount of charge 
transfer can be calculated from the change in the Fermi surface area including the SS, outer 
and inner Rashba bands using Luttinger’s theorem. 27 The charge transfer of Cr atoms after 
46 minutes Cr deposition is estimated to be q ~ 0.06 e~ per unit cell. The Rashba coupling 
parameter a can be extracted using a = A/cj?/2m*, where the m* is the effective mass of 
2DEG. From Fig. [2^1) , the maximum coupling parameter is 0.34 eV-A, which is comparable 
with those of Au(lll) surface^ 28 29 and 2DEG in semiconductors.^^ 

The absence of gap opening at the Dirac point and the emergence of Rashba-splitting 
bands suggest that the Cr atoms deposited on the surface interact with TI very differently 
from bulk doping. — In order to understand why surface deposition of Cr atoms does not lead 
to the gap opening as in the case of bulk Cr doping, we use DFT to simulate the deposition 
process of Cr atoms on Bi 2 Se 3 (llI) surface. At the beginning of the simulation, we put 
Cr atoms on different initial positions (including intercalation site, hollow site, bridge site, 
Bi-top site and Se-top site), and let them fully relax. It is found that Cr atoms on the Se-top 
site and bridge site will spontaneously relax to the Bi-top site. Therefore, we only have 
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Figure 3: (Color online) (a) A schematic diagram showing the deposition process of Cr atom 
to the surface of Bi 2 Se 3 . (b-f) Binding energies of Cr absorption on different sites: (b) hollow; 
(c,d) intercalation: (e) Bi-top; (f) Cr-dimer absorption on the hollow site. (b,c,e,f) are the 
top views and (d) is the side view. Black solid lines show the binding energies of absorption 
on different sites. The blue dash dot line stands for the chemical reaction path from climbing 
NEB calculation. The energy barrier is 0.907 eV. 

three energetically stable configurations as shown in Fig. [3j They are Cr absorptions on 
the hollow site [Fig. J3|b)], intercalation site [Fig. [3](c,d)], and Bi-top site [Fig. |3](e)], with 
binding energies (E^)^ of -2.378 eV, -2.631 eV and -2.200 eV respectively. Because E b of all 
optimal structures are negative, the energy will be released during the deposition process, 
and the structures with absorbed Cr atoms are stable. Although the configuration for Cr 
atom on the intercalation site has the lowest E^, this structure is practically difficult to 
be realized. According to the climbing nudged elastic band (NEB) calculation, 3 ^ a large 
energy barrier (about 1 eV) should be overcome when Cr atom intercalates into the atomic 
layers underneath. Therefore, the most feasible configuration observed experimentally is Cr 
absorption on the hollow site with the relatively lower E^. This is very different from the 
substitutional Bi site for bulk Cr doping. 15 l 

Here, we focus on the electronic structure of Cr-single absorption on the hollow site, and 
its density of state (DOS) is shown in Fig. |dj Due to C^ v symmetry around the absorption site, 
the crystal field splits Cr 3d orbitals into A\ (d z 2 ), E\ ( d XZl d yz ) and E 2 ( d x 2 _ y 2 , d xy ) states. 
Considering the Hunt’s rule, a large energy splitting between the spin-majority states and the 
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spin-minority states can be observed from our DFT calculations, and the five spin-majority 
states of Cr (Aj, e\ and e\ ) are fully occupied [see Fig. 4ja)]. Due to the similar in-plane 
orientation, a strong p-d hybridization between Se p orbitals and Cr d x 2_ y 2 and d xy orbitals 

_A* 

drives ^ states close to the Fermi level, and the spin-polarized electron density becomes 
more extended which also could be found on the neighboring Se atoms. Therefore, the total 
magnetic moment of this structure is 5 /ig. 

For magnetic properties of Cr absorption on the hollow site, we calculate the magnetic 
anisotropy energy (MAE) from DFT simulations. It is found the magnetic easy axis is along 
the in-plane direction with MAE of 2.04 meV/atom. Based on the perturbation theory G5EE3 
around the Fermi level, MAE is mainly determined by the difference of SOC matrix elements 
= (< \H so (n)\j)) along in-plane and out-of-plane direction, which can be expressed 
as — \Hij°(z)\ 2 ), where i and j stand for the occupied and un-occupied 

states respectively and Ay is their energy difference. Shown in Fig. |f][b), around the Fermi 
level, the states of Cr atom are mainly contributed by hybridized d orbitals with opposite 
spins, and the coupling between them results in an in-plane magnetic easy axis. 35 36 In TIs, 
the SS is protected by TRS and can only be gapped by out-of-plane ferromagnetism 8 or 
in-plane ferromagnetism combined with mirror symmetry breaking. 37 The absence of gap 
opening at the Dirac point in the ARPES data is in agreement with the in-plane magnetic 


easy axis. 19 

If absorbed Cr atoms form in-plane ferromagnetism without breaking other symmetry, 
the effective exchange held will shift the Dirac point of SS without opening a band gap. 
The shift magnitude (k®) can be estimated by kP = m%/vp, in which m% is the in-plane 
ferromagnetic coupling strength and vp is the Fermi velocity of BDSes (vp = 3.55 eV -A 38 ). 
To estimate the shift of Dirac point, in our view, the maximum magnetic coupling strength 
should have the same order of magnitude with that of Cr-dimer absorbed on Bi 2 Se 3 (lll) 
surface. Considering Cr-dimer absorption on the hollow site, calculated results are shown 
in Fig. [3|f): the distance between two absorbed Cr atoms is 2.87 A and E^ is -2.409 eV. 









Energy (eV) 

Figure 4: (Color online) (a) The total and orbital-projected density of states (DOS) of Cr 
atom for the single adsorption on the hollow site, (b) The a x component of DOS. The positive 
and negative stand for the spin-majority and spin-minority, respectively. The Fermi level is 
set to zero and marked by the red dot line. 
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This value is slightly lower than that of single Cr absorption on hollow site [see Fig. [3](b)], 
suggesting both structures can be observed in the experiment. The magnetic easy axis of Cr 
dimer is also along the in-plane direction with MAE of 3.32 meV/atom, and two Cr atoms 
prefer ferromagnetic coupling with coupling strength of 69.21 meV, and the average magnetic 
moment of each Cr atom is about 3.2 /ig. If the in-plane ferromagnetic coupling formed in Cr 
deposited Bi 2 Se 3 surface, using kP = rn^jv-p , we estimate the maximum momentum shift k^ 
is about 0.02 A” 1 , which is within our experimental uncertainty. We note that the absence of 
long range ferromagnetic order was reported 19 in previous studies of surface deposition of Fe 
on Bi 2 Se 3 surface. Here, we tend to believe that there is no long-range ferromagnetic order 
down to 15 K for surface deposition of Cr, which needs to be confirmed by direct magnetic 
measurements. 

In conclusion, our combined experimental and theoretical study show that, different from 
bulk Cr doped Bi 2 Se 3 , surface deposition of Cr atoms on Bi 2 Se 3 shows no gap opening at 
Dirac point, implying that no out-of-plane ferromagnetism is formed on surface down to 
our measurement temperature of 15 K. According to DFT calculations, the difference of 
electronic and magnetic properties originates from the different Cr configurations: surface 
deposition on the hollow site versus substitutional site bulk doping. Due to the coupling 
between Cr 3d orbitals around the Fermi level, the magnetic easy axis of deposited Cr atoms 
is along the in-plane direction. In addition, deposited Cr atoms induce electron doping to 
the sample and form Rashba-splitting 2DEG bands near the sample surface. These results 
indicate the interaction between magnetic dopants and TIs strongly depends on the intricate 
configuration of magnetic atoms, which has important implication for manipulation of SS in 
TIs. 
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Methods 

ARPES measurements were carried out at Beamline 12.0.1 of the Advanced Light Source 
with total energy resolution better than 25 meV. High quality bulk single crystal BDSes 
samples were cleaved in situ in UHV chamber with base pressure of 2.5xl0 -11 Torr before the 
measurements. The sample was kept at 15 K during Cr deposition and ARPES measurements. 
Cr was evaporated from a resistively heated tungsten basket, and the flux was calibrated by 
a quartz microbalance to be ~ 0.17 A per minute prior to experiment. Due to the different 
sticking coefficients between quartz and Bi 2 Se 3 (lll) surface, the actual amount of Cr atoms 
absorbed on Bi 2 Se 3 surface cannot be obtained directly. In order to estimate the amount 
of Cr absorbed on BDSes surface, we compare the charge transfer from the Bader charge 
analysis® 9 in the DFT calculation with the experimental charge transfer. The Bader charge 
analyis shows each Cr atom transfers 0.361 electron to B^Ses surface. Experimentlly, using 
Luttinger’s theorem, 27 the total charge transfer from Cr atoms calculated from the change 
of Fermi surface area after deposition for 46 minutes is 0.06 electron per unit cell. Thus we 
estimate 0.06/0.361 ~ 1/6 Cr atom was absorbed per unit cell after 46 minutes deposition, 
and the corresponding coverage is ~ 9% ML Cr. 

DFT calculations were carried out by using the projector augmented wave method®®] 
and the generalized gradient approximation with Perdew-Burke-Ernzerhof type functional, 42 
as implemented in the Vienna ab initio simulation package. 43 Plane wave basis set with a 
kinetic energy cutoff of 340 eV was used. For the geometrical and electronic structures, a 
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3x3x1 supercell slab model with the in-plane experimental lattice constants was constructed, 
which contained 20 atomic layers (4 QL) of Bi 2 Se 3 , absorbed Cr atoms, and a vacuum 
layer larger than 20 A along the z direction. During the structure optimization, 16 atomic 
layers of B^Ses on the bottom were fixed and the other atoms were fully relaxed until the 
residual force less than 0.01 eV /A. The Monkhorst-Pack k points were 2x2x1. SOC was 
included to calculate the electronic structure. 44 For the climbing nudged elastic band (NEB) 
calculations,^ 3 the 3x3x1 supercell slab model contained 5 atomic layers (1 QL) of Bi 2 Se 3 , 
absorbed Cr atoms and a vacuum layer with fully relaxation. 
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